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We present the design for an absorbing metamaterial element with near unity absorbance. Our 
structure consists of two metamaterial resonators that couple separately to electric and magnetic 
fields so as to absorb all incident radiation within a single unit cell layer. We fabricate, charac- 
terize, and analyze a metamaterial absorber with a slightly lower predicted absorbance of 96%. 
This achieves a simulated full width at half maximum (FWHM) absorbance of 4% thus making 
this material ideal for imaging purposes. Unlike conventional absorbers, our metamaterial consists 
solely of metallic elements. The underlying substrate can therefore be chosen independently of the 
substrate's absorptive qualities and optimized for other parameters of interest. We detail the de- 
sign and simulation process that led to our metamaterial, and our experiments demonstrate a peak 
absorbance greater than 88% at 11.5 GHz. 



The nascent field of electromagnetic metamaterials has 
produced exotic effects such as negative index of re- 
fraction [TJ [2], and devices such as an electromagnetic 
cloak [3]. The realization of such properties lies in the 
ability of metamaterials to create independent tailored 
electric [I] and magnetic [5J responses to incident radia- 
tion. Electromagnetic metamaterials are also geometri- 
cally scalable which translates into operability over a sig- 
nificant portion of the electromagnetic spectrum. To date 
metamaterials have been demonstrated in every tech- 
nologically relevant spectral range, from radio [DJ, mi- 
crowave |7], mm- Wave [8], THz [DJ, MIR [ID], NIR [IT], 
to the near optical [T2]. These designer electromagnetic 
materials are an ideal platform for the investigation of 
novel emergent physical phenomena while also holding 
great promise for future applications. 

As an effective medium [13] . metamaterials can be 
characterized by a complex electric permittivity e(w) = 
61 + ic2 and magnetic permeability jx{ui) = \x\ + ifj®. 
Much of the work in metamaterials has been focussed 
on the real part of e and /i to enable the creation of a 
negative refractive material. However, the oft-overlooked 
loss components of the optical constants (e2 and ^2) have 
much potential for the creation of exotic and useful ma- 
terials as well. For example, they can be manipulated 
to create a high absorber. By manipulating electric and 
magnetic resonances independently, it is possible to ab- 
sorb both the incident electric and magnetic field. Ad- 
ditionally, by matching e and /x, a metamaterial can be 
impedance-matched to free space, minimizing reflectiv- 
ity. In this paper, we show that metamaterials can be 
fashioned to create narrow-band perfect absorbers, with 
the potential to be used in devices such as in bolometers. 

A single unit cell of the absorber consisted of two dis- 




FIG. 1: Electric resonator (a) and cut wire (b). Dimension 
notations are listed in (a) and (b). The unit cell is shown 
in (c) with axes indicating the direction of propagation of a 
TEM wave. 



tinct metallic elements as shown in Fig. [T] (a) and (b). 
Electric coupling was supplied by the electric ring res- 
onator (ERR) , similar to the design presented by Padilla 
et al. [14], and is shown in Fig. [l](a). This element con- 
sisted of two standard split ring resonators connected by 
the inductive ring parallel to the split-wire. We used this 
design instead of a conventional split-wire design because 
of the limitations of straight wire media [IS], i.e. a split- 
wire design has limited tunability beyond the addition of 
more wires per unit cell to increase inductance |16j . The 
magnetic coupling required a more complicated arrange- 
ment, and thus in order to couple to the incident H- field, 
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FIG. 2: Simulated perfect metamaterial absorber. Re- 
flectance (green curve) and absorbance (red curve) are plotted 
from zero to 100% (left axis). The transmission is plotted on 
the right axis as the blue curve on a scale from zero to 5%. 



we needed flux created by circulating charges perpendic- 
ular to the propagation vector. We created this response 
by combining the center wire of the electric resonator 
with a cut wire (Fig. JTTb)) in a parallel plane separated 
by a substrate (Fig. fTfc)). This design is similar to the 
so-called "fishnet" and paired nanorod structures |lllll2j . 
in the sense that they derive a magnetic response from 
driving two antiparallel currents in conducting segments. 
We were then able to tune the magnetic response by 
changing the geometry of the cut wire and the separa- 
tion between the cut wire and electric resonator. By ma- 
nipulating the magnetic coupling without changing the 
geometry of the electric resonator, we able to decouple e 
and /i, and therefore individually tune each resonance. 

We performed computer simulations of an ideal - but 
realizable - metamaterial absorber using the commercial 
finite-difference time domain (FDTD) solver Microwave 
Studio by CST.[T7] The program simulated a single unit 
cell as shown in Fig. [I] (c) with appropriate bound- 
ary conditions, i.e. perfect electric (yz-plane) and per- 
fect magnetic (iz-plane). Waveguide ports on the other 
boundaries simulated a TEM plane wave propagating 
through the medium. Simulation produced the complex 
frequency dependent S-parameters, Su and S%i, where 
T(lu) — \S2i\ 2 and R(uj) = |Sii| 2 are the transmission 
and reflectance, respectively. From the S-parameter data 
we inverted the fresnel equations to extract the complex 
optical constants [TS]. We then examined the behavior of 
the surface current density, magnetic and electric fields 
at ujq to verify that we were coupling to the correct res- 
onant mode of each metamaterial element. By tuning 
ujq of the e and /i resonances (not shown), due to the 
two metamaterial elements, we achieved e = fi and thus 



an impedance near the free space value. While 100% 
asborbance is theoretically possible, it can only occur 
when the metamaterial layer is impedance-matched to 
free space such that the reflectance is zero.[T3] One may 
then minimize transmission such that the addition of 
multiple layers ensures T(u>) — > 0. The simulated meta- 
material had the dimensions, in millimeters, of: ai=4.2, 
£12=12, W=3.9, G=0.606, t=0.6, L=1.7, H=11.8, and the 
metamaterials elements were separated by 0.65 in the z 
direction. In Fig. [2] we show the simulation results for 
the metamaterial perfect absorber. Reflectance and ab- 
sorbance are plotted from zero to 100% (left axis) and 
the transmission is plotted on the right axis from zero 
to 5%. The reflectance is large ~ 97% near the bounds 
of the plot, 9GHz and 14GHz, but there is a minimum 
of 0.01% at uiq = 11.65GHz. The simulated transmis- 
sion also undergoes a minimum near ujq and yields a 
value of ~ 0.9%. Thus we achieve a best simulated ab- 
sorbance A{u>) = 1 — T(uj) — R(oj) slightly less than unity 
A(lu) = 99% with a FWHM of 4% compared to w . 

While our simulated absorber achieves narrow-band 
high absorbance, we were restricted by the minimum line 
widths available to us (250 /im) and other various fab- 
rication tolerances. These limitations were incorporated 
into the design process and we fabricated a metamate- 
rial which deviated slightly from the ideal absorber and 
had the dimensions, in millimeters, of: ai=4.2, a 2 =12, 
W=4, G=0.6, t=0.6, L=1.7, H=11.8. Each metalliza- 
tion was fabricated on a FR4 substrate with a thickness 
of 0.2mm. Metamaterials were fabricated using standard 
optical lithography. Photosensitized half-ounce copper- 
clad FR-4 circuit board formed the substrate material 
with an 17 /Ltm copper thickness. A mask was designed 
and printed in high resolution on a transparency. The 
photosensitized copper clad FR-4 was covered by the 
mask and exposed to ultraviolet light then developed in 
a sodium carbonate solution and etched in ferric chlo- 
ride. Finally, the board was placed in a stripper so- 
lution to remove the remaining photoresist. Both the 
ERR layer (Fig. [TJa)) and the cut-wire layer (Fig. JTTb)) 
were fabricated in such a fashion, see inset to Fig. pfrb). 
These boards where then sandwiched (using an adhesive 
with 0.06mm thickness) about another 0.2mm thick FR4 
blank substrate to obtain the correct spacing. The end 
results was that the metamaterials elements were sepa- 
rated by 0.72mm in the z direction. This design allowed 
us to couple to both magnetic and electric resonances 
without using out-of-plane elements typical with "wine- 
crate" designs, |2J thus greatly simplifying construction of 
the absorber. 

We experimentally verified the behavior of the ab- 
sorber by measuring the complex S-parameters of a large 
planar array of pixels (outer dimensions of 15 x 15 cm). 
We used an Agilent vector network analyzer that pro- 
duced microwaves in the range of 8-12 GHz. One mi- 
crowave horn (Rozendal and Associates) focused the GHz 
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FIG. 3: Simulated (red curves) and measured (blue curves) 
T(w) and R(w) for the microwave absorber, (a) displays the 
simulated and measured reflectance from zero to 100%. A 
schematic of the reflection experiment is shown as the inset 
to (a) where BS is the beam splitter. The simulated and 
measured transmission, shown in (b), are displayed from zero 
to 10%. Both reach a minimum of approximately 0.6% near 
11.5 GHz. Inset to (b) shows photographs of the individual 
components which make up the metamaterial absorber. The 
sample consists of the electric ring resonator (left) and the 
split wire (right) and is joined with an FR4 spacer of 0.72mm 
thickness. 



beam on the sample and another horn served as a detec- 
tor. Both horns coupled to linearly polarized light and 
had parallel polarization directions. To measure trans- 
mission, we set the horns in a normal incidence confocal 
configuration. We performed reflectance measurements 
at normal incidence. However, in order to eliminate large 
voltage standing wave ratio (VSWR) which would over- 
whelm the data, we used a beam splitter configuration, 
shown as the inset to Fig. [3] (a). Radiation first traveled 
from the transmitting horn through the beam splitter, 
reflected from the sample and returned to the beam split- 
ter where it was reflected at 45% to the receiving horn. 
We set the normalization and phase in the transmission 
configuration by simply removing the sample. For re- 
flection, we replaced the sample with a perfect reflector 
for normalized measurements. Transmission characteri- 
zation techniques utilized here have been described else- 
where in detail[20 . 

Computer simulations were performed for this meta- 
material structure and are presented for comparison to 
experimental data. The simulated R(uj) is shown as the 
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FIG. 4: Main panel shows the simulated (red) and measured 
(blue) absorbance curves. Significant broadening in the reflec- 
tion curve contributes to the deviation of the experimental ab- 
sorbance curve from simulation. The dashed gray absorbance 
curve is a Gaussian weighted average which approximates as- 
sembly error of the metamaterials absorber, described in the 
text. Inset shows the simulated angular dependence of the 
absorbance at at uimax- 



red curve in Fig. [3] (a). At low frequencies the reflectivity 
is high and yields values near 95%. A large feature oc- 
curs in R(uj) and has a minimum of ~ 3% near 11.5 GHz. 
Likewise the simulated transmissivity, shown as the red 
curve in Fig. |3](b), has 4% T(u) at low frequencies and 
then undergoes a minimum near 11.25 GHz. The simu- 
lated absorbance of our metamaterial unit cell thus has 
a maximum near where both R(lu) and T(lu) have their 
minima. We plot A(uj) as the red curve in Fig. [4] from 8 
to 12 GHz. The simulated A(u) peaks at 96% at 11.48 
GHz and has a FWHM of 4% with respect to its center 
frequency. 

By performing simulations for transmissivity and re- 
flectivity, we are able to account for the form of the ex- 
perimental data. The measured R(lu) and T(u>) are plot- 
ted as the blue curves in Fig. [3] (a) and (b) respectively. 
We find excellent agreement between simulated and mea- 
sured T(a>) over all frequencies characterized. The mea- 
sured reflectance is consistently 8% lower than the sim- 
ulated R(lo) at low frequencies. Both the simulated and 
experimental R(w) reach a minimum at approximately 
11.5 GHz, but experimentally the minimum is 11%, as 
opposed to the simulated value of 3%. There is also sig- 
nificant broadening in the experimental reflectivity curve 
with respect to simulation. This is most likely due to 
tolerances in the fabrication and assembly. 

The experimental absorbance calculated from T(lS) 
and R(u>) is shown as the blue curve in Fig. g The 
simulated and experimental curves have maximum ab- 
sorbance at the same frequency, w mQa; =11.5 GHz. How- 
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FIG. 5: Simulated absorbance with increasing metamaterial 
layers, red curve is a guide to the eye. Insets show the simu- 
lated losses at resonance. The Ohmic loss (surface) is shown 
in the left inset panel and the dielectric (volume) loss is on 
the right. 



ever, the simulated curve reaches a maximum of 96% and 
the experimental curve only reaches 88%. The differ- 
ence between the simulated and measured curves is due 
primarily to deviations in R{ui). Better fabrication and 
assembly tolerances will permit absorbances and band- 
widths similar to that shown as the red curve in Fig. 
[4] We performed computer simulations to verify the ori- 
gin of the broadening of A(u>). From a theoretical view- 
point it is expected that the frequency location of the 
absorbance peak depends significantly upon the separa- 
tion between the electric ring resonator and the cut-wire, 
since this distance determines the frequency location of 
the magnetic resonance. Computer simulations verified 
this conjecture. We simulated various spacings centered 
about the designed value of 0.72mm, and computed a 
Gaussian weighted average A(u>) with a standard devi- 
ation of a =20 fjxa. The absorbance calculated in this 
manner agrees very well with the experimental curve and 
is plotted as the grey dashed curve in the main panel of 
Fig. |4j Thus assembly errors in the metamaterial spac- 
ing of only 5% can easily account for the disagreement 
between the simulation and experiment shown in Fig. [4] 
as the red and blue curves, respectively. 

An investigation into the effect of adding multiple 
metamaterial absorbing layers is shown in the main panel 
of Fig. [5j As can be observed the absorbance rises 
sharply with additional layers and is asymptotic to unity, 
within computational noise. Two layers of the metama- 
terial absorber achieve a value of ^99.9972%. At this 
thickness the entire metamaterial is only ~ Ao/2 at res- 
onance. For the structure shown in Fig. [lj A(a>o) can 
never reach unity since the metamaterial is not perfectly 
matched to free space. However, it is important to note 



that although an impedance matched structure may yield 
zero R(wo)j this does not necessarily mean the trans- 
mission will be small. In this study we are interested 
in achieving R(w o )=T(u;o) = in a single unit cell in the 
propagation direction. Thus our metamaterial structure 
was optimized to maximize the absorbance with the re- 
striction of minimizing the thickness. If this constraint is 
relaxed, however, an impedance matched condition can 
be obtained and with multiple layers a perfect absorbance 
can be achieved. 

We further investigated the origin of the loss in our 
metamaterial structure, as indicated in the inset to Fig. 
[5] The Ohmic "surface" loss originates mainly from the 
center conducting region of the electric ring resonator, 
and dielectric losses primarily occur in-between the two 
metamaterial elements where the electric field is large. 
The major component responsible for the absorbance in 
our metamaterial structure is the dielectric loss and simu- 
lation indicates this is an order of magnitude greater than 
Ohmic losses. This is consistent with previous studies of 
frequency selective surfaces (FSSs) where it was found 
that metallic absorption was relatively insignificant in 
comparison to dielectric losses. [2"T] 

We now discuss the potential application of a perfect 
metamaterial absorber as a bolometer. An ideal bolo- 
metric detector is one in which all photons that fall upon 
its surface are absorbed, converted to heat, and thus 
sensed. 22J A prominent feature of commercially avail- 
able bolometers is their extreme bandwidth making them 
ideal candidates as detectors, especially in the far in- 
frared frequency regime where other efficient detection 
methods are largely unavailable. However, bolometers 
which achieve a narrowband response also have signifi- 
cant application as focal plane array detectors for imag- 
ing. The metamaterial design presented here achieves 
simulated absorbtivities near unity making them ideal 
candidates for bolometric pixel elements. The natu- 
ral narrowband metamaterial resonance is a salient fea- 
ture for focal plane array detectors as it is naturally 
apodizing and operates at room temperature. Since 
the elements which constitute our bolometer are sub- 
wavelength, metamaterials can inherently image at the 
diffraction limit. Further, here we have used only a sin- 
gle unit cell in the propagation direction, (with a thick- 
ness significantly smaller than the wavelength Ao /35) , yet 
achieved an experimental absorbance of 88%. By adding 
multiple layers one can achieve narrowband absorptivi- 
ties of unity. Although the design is planar and A(u) 
should fall off relatively rapidly, we investigated the an- 
gular dependence (inset to Fig. [4} and found the meta- 
material still achieved 50% absorbance at an incident 
full-angle of 16°. Finally, the scalability of metamateri- 
als permits the usage of these perfect absorbers at other 
wavelengths of interest such as the possibility of room 
temperature high resolution imaging at mm- Wave and 
THz frequencies. This is particularly interesting owing to 
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the many possible applications. We have designed simi- 
lar metamaterials (not shown) as in Fig. Illthat achieved 
comparable A{uS) to that shown in Fig. 15] operating at 
94 GHz and 1 THz. 

While the creation of a metamaterial absorber is novel, 
the technology to create highly absorptive materials in 
the microwave frequency range is well-established. The 
material "chiroshicld" is capable of very high absorption, 
and reduces backscatter by 15 to 25 dB [23] . Chiroshield, 
however, has geometries on the order of one wavelength, 
which is 35 times larger than our design at a given center 
frequency. Some ferrites show reflection loss of -30 dB 
[24] . These ferrites, however, lack the inherent and pre- 
cise tunability of metamaterials. Other work has been 
done on the application of frequency-selective materials 
to bolometric devices. These bolometric designs, how- 
ever, require cryogenic temperatures to operate and only 
achieve absorbances of about 50% [25 .. At optical wave- 
lengths, metal colloides are known to have large absorp- 
tive properties, [55] due to the geometric specific surface 
modes [57]. This suggests and interesting possible appli- 
cation of artificial plasmonic wire media^Q shaped into 
particular geometries to function as perfect absorbers [16] 
operating at microwave frequencies. 

In conclusion, we have demonstrated that metamate- 
rials can be highly absorptive over a narrow frequency 
range. This stems from the ability to design metamate- 
rial elements which can individually absorb the electric 
and magnetic components of an incident electromagnetic 
wave. In contrast previous experimental results at THz 
frequencies using a single type of metamaterial element 
(electric) yielded values of only A(w)=~ 20%. [28] The 
resonant frequency of the metamaterial elements may 
be tuned throughout some range of frequencies [35J [3D] 
thus enabling hyperspectral imaging. A further benefit 
afforded by metamaterials is the ability to construct a 
single unit cell with /j,(lu) = e(uj) over an extended fre- 
quency range. Thus this unit cell can achieve zero re- 
flectance since it can have an impedance equal to the 
free space value Z = \/ fi/e = 1. This is similar to the 
well known theoretical construct, the Perfectly Matched 
Layer (PML) Absorbing Boundary Condition (ABC).[3T] 
which splits waves incident upon a boundary into electric 
and magnetic components to obtain near perfect absorp- 
tion. The PML, however, requires gain and additionally 
is on the order of Ao/2 in thickness |32j. 

The design presented here can still be improved. By 
incorporating a substrate with a highly consistent dielec- 
tric constant, we will be able to optimize the design at 
the correct resonant frequency. Similarly, tighter fabri- 
cation tolerances would allow us to bring the impedance 
closer to unity and therefore closer to 100% absorption. 
More progress could also be made on the design itself. 
Currently, the absorber is polarization sensitive, which 
is not ideal for some applications. Additionally, the as- 
sembly is complicated due the restrictions placed on the 



direction of the incident wave to couple to the magnetic 
resonator. This resulted in variations between unit cells 
and less distinct absorbance peak. 

NIL and WJP acknowledge support from the Los 
Alamos National Laboratory LDRD program, DRS ac- 
knowledges support MDA Contract W9113M-07-C-0078, 
and DRS, SS, and JJM acknowledge support from a Mul- 
tiple University Research Initiative sponsored by the Air 
Force Office of Scientific Research Contract No. FA9550- 
06-1-0279. 



[4 
[5 

[6 

[7 
[8 

[9 

[10 
[11 

[12 
[13 
[14 

[15 
[16 
[17 
[18 
[19 

[20 
[21 



V. G. Vesalago, Sov. Phys. Usp. 10, 509 (1968). 

R. A. Shelby, D. Smith, and S. Schultz, Science 292, 77 

(2001). 

D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. 
Pendry, A. F. Starr, and D. R. Smith, Science 314, 977 

(2006) . 

J. B. Pendry, A. J. Holden, W. J. Stewart, and I. Youngs, 
Phys. Rev. Lett. 76, 4773 (1996). 

J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. 
Stewart, IEEE Trans. Microwave Theory Tech. 47, 2075 
(1999). 

M. C. K. Wiltshire, J. B. Pendry, I. R. Young, D. J. 
Larkman, D. J. Gilderdale, and J. V. Hajnal, Science 
291, 849 (2001). 

D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat- 
Nasser, and S. Schultz, Phys. Rev. Lett. 84, 4184 (2000). 
M. Gokkavas, K. Guven, I. Bulu, K. Aydin, R. S. Penciu, 
M. Kafesaki, C. M. Soukoulis, and E. Ozbay, Phys. Rev. 
B 73, 193103 (2006). 

T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. 
Smith, J. B. Pendry, D. N. Basov, and X. Zhang, Science 
303, 1494 (2004). 

S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, 

and C. M. Soukoulis, Science 306, 1351 (2004). 

S. Zhang, W. Fan, N. C. Panoiu, K. J. Malloy, R. M. 

Osgood, and S. R. J. Brueck, Phys. Rev. Lett. 95, 137404 

(2005). 

G. Dolling, M. Wegener, C. M. Soukoulis, and S. Linden, 
Opt. Lett. 32, 53 (2007). 

D. R. Smith and J. B. Pendry, J. Opt. Soc. Am. B 23, 
391 (2006). 

W. J. Padilla, M. T. Aronsson, C. Highstrete, M. Lee, 
A. J. Taylor, and R. D. Averitt, Phys. Rev. B 75, 041102 

(2007) . 

D. Schurig, J. J. Mock, and D. R. Smith, Appl. Phys. 
Lett. 88, 041109 (2006). 

D. R. Smith, D. C. Vier, W. Padilla, S. C. Nemat-Nasser, 
and S. Schultz, App. Phys. Lett. 75, 1425 (1999). 
CST Studio Suit e 2008, CST of America. 
http://www.cst.com (2008). 

D. R. Smith, S. Schultz, P. Markos, and C. M. Soukoulis, 
Phys. Rev. B 65, 195104 (2002). 

S. G. Johnson, P. Bienstman, M. Skorobogatiy, 
M. Ibanescu, E. Lidorikis, and J. D. Joannopoulos, Phys. 
Rev. E 66, 066608 (2002). 

A. F. Starr, P. M. Rye, D. R. Smith, and S. Nemat- 
Nasser, Phys. Rev. B 70, 113102 (2004). 
J. E. Raynolds, B. A. Munk, J. B. Pryor, and R. J. 



G 

Marhefka, J. Appl. Phys. 93, 5346 (2003). 
[22] P. L. Richards, J. App. Phys. 76, 1 (1994). 
[23] I. V. Lindell, A. H. Sihvola, S. A. Tretykov, and A. J. Vi- 

itanen, Electromagnetic waves in Chiral and Bi-Isotropic 

Media (Artech House, 1994). 
[24] P. Singh, V. K. Babbar, A. Razdan, S. L. Srivastava, and 

T. C. Goel, Materials Science and Engineering B 78, 70 

(2000). 

[25] T. A. Perera and et al., App. Opt. 45, 7643 (2006). 

[26] W. Lamb, D. M. Wood, and N. W. Ashcroft, Phys. Rev. 
B 21, 2248 (1980). 

[27] J. B. Pendry, J. Mod. Opt. 41, 209 (1994). 

[28] J. F. O'Hara, E. Smirnova, H.-T. Chen, A. J. Taylor, 
R. D. Averitt, M. L. C. Highstrete, and W. J. Padilla, J. 
Nanoelectron. Optoclcctron. 2, 90 (2007). 

[29] W. J. Padilla, D. R. Smith, and D. N. Basov., J. Opt. 
Soc. Am. B 23, 404 (2006). 

[30] H.-T. Chen, J. F. O'Hara, A. K. Azad, A. J. Taylor, 
R. D. Averitt, D. Shrekenhamer, and W. J. Padilla, Na- 
ture Photonics - To appear (2008). 

[31] J.-P. Berenger, J. Comp. Phys. 114, 185 (1994). 

[32] A. Taflove and S. C. Hagness, Computational Electro- 
dynamics: The Finite- Difference Time-Domain Method, 
Third Edition (Artech House, 2000). 



